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SFCTION  1 

SUMMARY 


Presented  herein  tire  the  results  of  ti  four-month  technical  program  to  develop 
and  demonstrate  the  application  of  an  inverse  perturbation  design  procedure, 
originally  developed  at  United  Technologies  Research  Center,  to  curved  structures. 
Briefly,  the  inverse  perturbation  technique  relies  on  information  concerning  the 
existing  vibration  modes  of  a structure  to  calculate  the  necessary  changes  to  ob- 
tain new,  prescribed  vibration  characteristics  for  selected  modes  of  the  structure, 
and  its  application  to  curved  structures  necessitates  the  use  of  finite-element 
computer  analyses  to  obtain  the  necessary  accuracy.  Consequently,  the  study- 
involved  formulation  of  the  equations  and  computer  programs  necessary  to  process 
NASTRAN  output  data  on  vibration  mode  shapes,  and  preparation  of  the  UTRC  inverse 
perturbation  design  program  to  accept  the  processed  output,  and  compute:  a)  the 
perturbation  functions  associated  with  changes  in  frequency  and  mode  shape  admixture 
( admixture  between  mode  shapes  implies  that  the  mode  shapes  of  an  altered  structure 
may  be  expressed  as  linear  combinations  of  the  mode  shapes  of  an  unaltered  structure), 
and  b)  the  perturbations  that  result  from  design  changes,  as  well  as  perform  checks 
on  the  stiffness  orthogonality  and  Rayleigh  quotient. 


Supported  under  Contract  F33615-77-C-2092 , sponsored  by  the  Air  Force  Aero 
Fropulsion  Laboratory,  Air  Force  Systems  Command,  the  program  successfully  demon- 
strated concept  feasibility  with  the  redesign  of  a flat  plate  and  a arc  of  a 
cylindrical  shell,  both  mounted  as  cantilever  beams.  Presented- in  the  following 
sections  are:  an  introduction  which  briefly  reviews  the  motivation  for  the  inves- 
tigations (Section  II ) ; an  overview  of  the  concept  of  inverse  perturbation 
(Section  III);  a description  of  the  computer  programs  written  to  implement  the 
concept  (Section  IV);  a presentation  of  the  computations  and  experiments  performed, 
together  with  the  results  obtained  for  the  two  cantilevered  objects  (Section  V); 
and  a discussion  of  potential  benefits  to  be  derived  from  continued  development  of 
the  technique  (Section  Vi).  Included  as  an  Appendix  is  a detailed  description  of 
previous  UTRC  work  on  redesigning  structural  vibration  modes  by  inverse  perturba- 
tion, subject  to  minimal  change  theory,  which  set  the  stage  for  the  investigations 
reported  herein. 
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SECTION  II 
INTRODUCTION 


The  design  of  structures  to  withstand  the  vibratory  excitation  experienced 
inside  a Jet  engine  is  clearly  a multifaceted  problem.  Nut  only  must  resonar.t 
response  be  considered,  but  also  the  possibility  of  aerodynamical ly  excited  vibra- 
tions such  as  flutter.  Furthermore,  perf  rmanee  requirements,  such  as  the  ability 
of  a fan  to  compress  air  or  of  a turbine  to  extract  energy  from  an  air  flow,  and 
practical  issues,  such  as  weight  and  resistance  to  impact  damage,  also  enter  into 
the  choice  f a good  design.  The  scope  of  the  problem  makes  a truly  general-purpose 
computerized  design  procedure  unlikely;  however,  situations  do  arise  where  answers 
to  a portion  of  the  total  problem  can  be  useful  in  formulating  the  proper  design 
approach. 

Among  the  most  enigmatic  facets  of  a structure  are  its  vibration  modes.  These 
are  functions  of  the  structure  that  can  be  used  to  characterize  its  response  to 
forced  vibratory  excitation  (Ref.  l),  to  transients,  and  even  to  static  forces 
IRef.  2).  So  central  are  these  functions  in  characterizing  vibratory  response, 
that  in  recent  years  numerous  large-scale  computer  programs  have  been  developed  to 
analyze  structures,  by  finite-element  methods,  to  determine  these  modes.  From  the 
experimental  point  of  view,  holographic  vibration  analysis  has  greatly  facilitated 
the  determination  of  vibration  modes  on  the  visible  surfaces  of  structures  (Ref.  1). 

Whereas  finite-element  and  holographic  methods,  in  addition  to  numerous  other 
experimental  and  theoretical  techniques,  make  it  possible  to  analyze  a given  struc- 
ture for  its  vibration  modes,  the  question  of  how  that  structure  should  be  changed 
to  effect  a specified  change  in  its  vibration  modes  lies  relatively  unapproached. 

Work  at  the  United  Technologies  Research  Center  in  recent  years  (Refs.  3,  k and  5) 
has  suggested  a potentially  powerful  tool,  based  upon  perturbation  equations,  for 
just  such  a task.  It  is  possible,  of  course,  to  use  perturbation  theory  to  gain  an 
approximate  answer  to  the  question  of  how  the  vibration  modes  of  a structure  will 
change  as  a result  of  small  changes  in  the  structure;  this  has  been  known  since 
the  19th  century.  The  interesting  problem  is,  however:  how  does  one  invert  the 
usual  perturbation  analysis  to  answer  the  question  posed  at  the  beginning  of  this 
paragraph?  The  most  recent  work  in  this  area,  at  these  laboratories,  has  established 
that  by  minimizing  the  change  to  a structure  necessary  to  accomplish  a set  of 
changes  to  its  vibration  modes,  we  may  define  a unique  set  of  functions  which  in 
turn  generate  a unique  change  to  the  structure.  These  functions  may  be  called  per- 
turbation functions,  and  they  are  formed  from  pairs  of  vibration-mode  functions  and 
their  spatial  derivatives.  A detailed  description  of  that  design  technique,  which 
we  call  inverse  perturbation,  together  with  its  application  to  the  redesign  of  a 
cantilever  plate,  based  upon  vibration  mode  data  obtained  from  hologram  interferometry 
is -included  as  Appendix  I to  the  present  report. 
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Presented  in  the  main  body  of  the  current  report  is  a description  of  the 
program  of  research  carried  out  under  Air  Force  Contract  F33615-77-C-0292  to 
extend  this  inverse  perturbation  method  of  design,  developed  for  holographically 
determined  vibration  modes,  so  that  it  cun  make  use  of  the  NACTRAN  finite-element 
program  for  vibration  analysis.  This  has  resulted  in  a design  program  which,  in 
its  present  form,  can  be  used  in  conjunction  with  NASTRAN  to  alter  the  thickness  of 
any  shell-like  structure  in  a vuy  that  will  generate,  approximately,  some  specified 
change  in  the  vibration  modes  of  a structure.  At  present,  changes  in  lateral 
dimensions  of  the  shell,  or  in  the  boundary  conditions  of  the  shell  cannot  be  dealt 
with,  but  future  programs  can  be  conceived  to  make  this  possible.  The  current  pro- 
gram is  likely  to  require  a number  of  iterative  designs,  each  involving  a NASTKAN 
analysis  to  finally  close  upon  a structure  whose  actual  vibrations  are  the  ones 
desired;  however,  all  indications  are  that  the  procedure  will  converge. 


SECTION  III 


THE  CONCEPT  OF  INVERSE  PERTURBATION 


The  details  of  the  derivations  involved  in  inverse  perturbation  can  be  found 
in  Rel's.  1*  and  5,  as  ■well  as  Appendix  I to  the  present  report.  Presented  here  are 
the  general  results  and  the  specific  equations  that  relate  to  this  program.  First, 
the  design  variable  is  taken  to  be  the  percentage  change  in  thickness  of  the  plate 
or  shell,  Ah(x,y )/h(x,y ) , where  h(x,y)  is  the  thickness  of  the  structure,  and  x and 
y are  spatial  coordinates  in  its  surface.  This  design  variable  is  expanded  as  a 
series  of  functions,  6(x,y),  called  perturbation  functions,  which  are  functions  of 
the  coordinates,  x and  y.  The  coefficients  of  the  series  are  determined,  in  part, 
by  the  set  of  changes  that  are  desired,  (A),  which  can  be  written  as  a row  matrix. 
The  (A)'s  are  related  to  percentage  changes  in  frequency,  and  to  coefficients 
of  admixture  between  mode  shapes.  (As  noted  in  the  Summary,  admixture  between  mode 
shapes  implies  that  the  mode  shapes  of  an  altered  structure  may  be  expressed  as 
linear  combinations  of  the  mode  shapes  of  an  unaltered  structure.)  The  final  set 
of  factors  in  the  equation  that  determines  the  percentage  changes  in  thickness  is 
the  inverse  of  what  we  may  call  a perturbation  matrix,  [Rj.  Thus  the  primary  equa- 
tion of  inverse  perturbation  is 


Ah( x,y )/h( x ,y ) = ( A ) [ B ] i(0(x,y))1.  (1  ) 

In  Eq.  (l),  the  necessary  set  of  perturbation  functions  have  been  arranged  as 
a row  matrix,  (0(x,y)),  and  the  superscript  T denotes  the  transposed  matrix;  i.e., 
a column  matrix.  In  order  to  define  the  matrix  [ B ] , it  is  necessary  to  specify  the 
set  of  changes  (a)  and  the  set  of  perturbation  functions  (0).  (Note:  the  functional 
dependence  of  0 on  x and  y is  implied  throughout.)  The  perturbation  functions  are 
formed  from  products  of  parameters  obtained  from  pairs  of  vibration  modes,  or  from 
squares  and  products  of  parameters  obtained  from  a single  mode.  Thus,  we  may  identify 
perturbation  functions  by  a double  subscript  0 or  0nn,  depending  upon  whether 
the  perturbation  function  is  formed  from  two  vibration  mode  shapes  <t>  ( x , y)  and 
♦jl(x*y),  or  from  <t’n(x,y)  alone.  Similarly,  the  change  parameters  (A)  may  be 
characterized  by  double  subscripts;  i ^ being  related  to  the  amount  of  mode  shape 
't’k  that  is  found  added  to  $u  as  a result  of  the  structural  change,  and  A,in  being 
related  to  the  change  in  frequency  of  the  nth  mode.  The  specification  of  any  change 
parameter  requires  the  inclusion  of  the  corresponding  perturbation  function  in 
Eq.  (l).  In  this  way , the  number  of  degrees  of  freedom  in  changing  the  structure 
(i.e.,  the  perturbation  functions)  matches  the  number  of  constraints  (i.e.,  the 
change  parameters ) . 
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With  this  in  mind,  let  us  define  the  elements  of  the  perturbation  matrix,  b . 
'l’iie  double  subscripts  of  tiie  change  variables,  and  the  perturbation  functions,  may 
be  ordered  in  some  arbitrary,  but  mutually  consistent  way.  Let  the  double  subscript 
'uk'  correspond  to  the  single  index  'p',  and  some  other  subscript  pair  ' lm'  corre- 
spond to  ' q'.  Then,  the  element  b is  defined  by 

bpq  = "ep(x,y)0q(x,y)  d*dy  . (2) 

where  the  integration  is  carried  out  over  the  surface  of  the  structure.  In  this 
particular  program,  however,  the  structure  lias  been  divided  into  discrete  elements 
so  tliat  the  integrations  will  extend  only  over  the  elements,  and  the  contributions 
of  the  elements  summed  to  obtain  tiie  coefficient  b 

pq 

In  order  to  define  the  perturbation  functions,  we  must  first  define  tiie  mode 
functions  according  to  the  form  they  take  in  the  NASTRAN  program.  Each  element 
of  the  shell-like  structures  considered  in  this  program  is  a triangular  platelet  of 
uniform  thickness  which  may  bend  out  of  its  plane  to  any  shape  describable  by  a 
cubic  polynomial  (minus  tiie  x“"y  term),  and  may  deform  and  rotate  homogeneously  within 
its  plane.  In  terms  of  a coordinate  system  local  to  each  element,  we  may  define 
tiie  mode  function,  vectorially,  as 


= i +1<t  + k $ 

A nx  J ny  nz ! 


where  i,  J and  k are  unit  vectors  in  the  x,  y and  z directions,  witli  z being  out 
of  the  plane  of  the  element,  and  the  origin  located  at  one  of  the  corners  of  the 
element.  The  components  of  4>n  are,  in  their  polynomial  forms: 


♦nx  = anl  + an2x  + an3y’ 


♦ny  = bnl  + bn2x  + 


♦nz  " cnl  + cn2x  + cn?  + cnUx  + en5^  + W + 

cnTx3  * cn8xy2  + W3- 


The  perturbation  functions,  which  may  now  be  defined,  are  most  easily  expressed 
as  the  sum  of  three  terms:  l)  an  inertial  term,  O ; 2)  a term  related  to  bending 
strain,  0 ; and  3)  a term  related  to  membrane  strain,  l\,m-  Thus, 


0 + 0 , + 0 
mp  sbp  snip  ’ 
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where  the  subscript  ’ p * denotes  the  mode  pair  ' nK * . The  perturbation  function, 
it  should  be  noted,  is  a scalar.  The  first  term  of  Eq  (7)  Is 


°mp  = -(tVh/t0kMk)  *n  ' \ ’ 


where  u>  and  10.  are  the  natural  frequencies,  in  radians,  of  modes  n and  k,  p is 
n k 

volume  density,  and  M is  the  model  mass  of  the  kth  mode,  defined  by 

K 

Mk  = // Ph  \ ' *k  dx  ^v-  ,t 


To  define  the  second  two  terms,  let  ua  first  define  the  matrix  G a 

'WJ 

1 v 0 

= t Kh/u>niokMk(l-V  )]  v 1 0 , 

0 0 (l-v)/2 


where  E is  Young’s  modulus  and  v is  Poisson's  ratio.  We  may  now  define  G^us 

c.  = G h2/ii 

%D  <\jS 


where , 


°sbn  - (♦„)  Sb 


(>1  ) = Uxx,  iyy  , ) , and 

kz  kz  kz  kz 


(*  ) = <♦**,  2<fx,v ) , 

nz  nz  nz  nz 


with  the  superscripts  denot.  inf,  partial  di  fferent.  iat.  ion  with  respect  to  the  super- 
script variables.  Finally, 


0 = (<t>  ) G (<f  ), 

smp  nxy  a.s  kxy 


where1 , 


(■I''  ) = (>!>x  , $y  , + <’•'),  and 

nxy ' nx  ny  ny  nx 


(Vv}  = Kx'  *ky ' *ky  v ^kx}- 

Now  let  ur.  define  the  change  parameters,  (A),  in  terms,  of  the  percentage 
changes  in  frequency  and  coefficients  of  admixture: 


(16) 


Ann  - [(1  + Awn/“n)2  - 1]  , and 

Ank  = cnk(a)n  ~ “k>/“n“k» 

where  is  the  admixture  coefficient.  The  new  mode  shape,  $ , is  expressible  in 

terms  of  the  old  mode  shapes  by  the  series  — 

*n  = + ZCnk  (18) 

where  the  term  for  n = k is  omitted  from  the  summation. 

The  format  of  the  NASTRAN  analysis  requires  some  degree  of  approximation  in 
implementing  Eq.  (l).  First  of  all,  whatever  new  design  is  generated  it  should 
have  plate  elements  of  constant  thickness,  or  else  the  analysis  of  the  new  struc- 
ture becomes  very  difficult.  This  means  that  the  values  of  the  perturbation  func- 
tions in  Eq.  (l)  should  be  made  constant  for  each  element  by  taking  the  average 
value  over  the  element.  Another  difficulty  arises  in  the  computation  of  the 
matrix  coefficients,  b , via  Eq.  (2).  If  the  integrations  over  each  element  are 
carried  out  as  indicated,  with  products  of  the  polynomials  of  Eqs . (4-6),  then 
the  resulting  integrands  may  be  as  high  as  twelfth  order  polynomials.  To  avoid  the 
awkwardness  of  having  to  deal  with  such  large  polynomials,  we  decided  to  approxi- 
mate the  perturbation  functions  within  the  elements  by  their  average  values  therein. 
This  made  evaluation  of  the  matrix  coefficients,  via  Eq.  (2),  simpler,  and  simul- 
taneously provided  constant  changes  in  thickness  for  the  plate  elements. 
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SECTION  IV 


COMPUTER  PROGRAM  DESCRIPTION:' 


A number  of  computer  programs  were  required  in  order  to  implement  the  design 
concept  of  the  subject  contract . As  described  in  the  following  paragraphs,  they 
included  the  primary  routines  for:  l)  accepting  the  NASTRAN  data,  processing  it 
and  forming  coefficient  data  files  (POSTPR);  C)  accepting  the  data  files  and  com- 
puting the  associated  perturbation  functions  ( NAGTY ) ; and  3)  accepting  the  per- 
turbation functions  and  performing  the  redesign  (CHANGNAGT),  as  well  as  various 
ancillary  routines  for  performing  checks  (TESTNASTY,  CP  and  DI8PPV ) , and  subroutines 
for  simplifying  some  of  the  programming. 


POSTPR  Program 

This  program  accepts,  as  input,  the  displacements  and  rotations  of  the  grid 
points  of  ttie  finite  elements  in  a NASTRAN  vibration-mode  analysis,  in  what  is 
referred  to  as  the  global  coordinate  system.  These  displacements  and  rotations  are 
converted  to  sets  of  coordinates  that  are  local  to  each  element.  The  displacements 
and  rotations  are  then  used,  together  with  Eqs.  ( U— 6 ) , to  solve  for  the  coefficients, 
that  determine  the  internal  deformations  of  the  elements;  i.e.,  the  u's , b’s,  and 
c's  of  Eqs.  ( U— 6 ) . This  is  done  for  each  NAGTRAN  vibration  mode,  and  the  resulting 
arrays  of  coefficients  are  stored  in  data  files  for  use  in  subsequent  programs. 


NAGTY  Program 


7h is  program  forms,  perturbation  functions  from  Eqs..  (i'-l  ),  and  tak< their 
average  value  over  each  triangle.  Two  subroutines.  (POLVOL  and  TKIAHG)  were  written 
to  facilitate  this  computation.  POLVOL  finds  the  value  of  the  integral  of  t lie  pro- 
duct of  two  polynomials  over  a triangular  region  in  the  x,  plane,  while  TK1ANG, 
a subroutine  to  POLVOL,  integrates  any  term  of  the  two-dimensional  polynomial  over 
the  triangular  region.  NAGTY  itself  is  a subroutine  to  the  main  inverse  perturba- 
tion progmm  (CIIANGNAST) . 


CHANGNAGT  Program 


This  is  the  primary  inverse  perturbation  design  program,  and  it 
written  no  that  it  may  use  perturbation  functions  generated  either  by 
therefore  NASTRAN)  or  by  PKRTNAGT;  an  updated  version  of  the  original 
generates,  perturbation  functions  from  experimental  holographic  data. 


has  been 
NAGTY  (and 
program  that 
CHANGNAGT 


8 


accepts  the  various  material  parameters  (Young's  modulus,  Poisson's  ratio,  density, 
and  the  vibration  mode  frequencies),  as  well  as  the  desired  change  parameters,  i.e., 
the  values.  It  thereupon  computes  the  necessary  perturbation  functions,  forms 
the  perturbation  matrix,  and  generates  the  new  design  via  Eq.  (l),  which  it  can  make 
available  on  punched  cards.  It  also  prints  the  modal  masses,  the  perturbation 
matrix,  the  determinant  of  the  perturbation  matrix,  the  coefficients  of  the  series  of 
the  perturbation  functions  that  comprise  the  new  design,  and  the  root -me an -square 
change  in  the  element  thicknesses.  Finally,  via  a subroutine  called  NASTDL,  it 
computes  the  various  perturbations  that  may  be  expected,  even  for  the  parameters 
that  are  unconstrained,  by  a formula  that  is  more  accurr.tc  than  the  perturbation 
formula  used  for  the  inversion  process  (see  Appendix^. 


TESTNASTY  Program 

This  program  will  print  the  values  of  the  perturbation  functions  formed  by 
NASTY,  and  compute  the  Rayleigh  quotients  (the  ratio  of  maximum  potential  energy  to 
maximum  kinetic  energy  in  the  vibration  cycle)  and  the  check  for  orthogonality  of 
the  stiffness  functions  (see  Appendix  ).  It  will  also  compare  NASTRAN  computations 
of  the  vibration  modes  with  holographic  data  when  the  vibration  mode  is  a scalar 
function.  It  does  this  by  computing  the  mode-function  value  of  the  NASTRAN  solution 
at  the  center  of  each  section  of  the  structure  for  which  the  holographic  data  has 
been  digitized. 


C2  Program 

This  program  takes  NASTRAN  solutions  for  the  vibration  modes  of  an  original 
and  modified  structure,  and  computes  the  admixture  coex’ficients  that  best  describe 
the  new  modes  as  a series  of  the  old  modes. 


DISP2V  Program 

This  program  computes  vectorial  displacements,  in  one  plane,  of  a plate  or 
shell  from  data  provided  by  two  holographic  interferograms  of  a vibration  mode. 
(The  third  component  is  assumed  to  be  zero.)  By  selecting  the  centers  of  the  sec- 
tions of  the  interferograms  to  be  digitized  so  that  they  correspond  to  the  NASTRAN 
grid  points,  direct  comparison  can  be  made  to  the  NASTRAN  analysis. 


SECTION  V 


ANALYTICAL-EXPER1 MENTAL  INVESTIGATIONS 


NASTRAN  Modeling  of  the  Vibration  Modes 

The  first  computations  performed  were  the  NASTRAN  analyses  of  both  a uniform 
plate  and  a uniform  shell  segment.  The  plate  was  15.2b  cm  long  by  IS. 7 em  wide  by 
0.317  cm  thick,  and  the  shell  segment  was  a 1*5°  arc  of  a cylindrical  shell,  7.62  cm 
long,  with  an  outer  radius  of  7-t>2  cm  and  a thickness  of  0.317  cm.  Schematic 
illustrations  of  the  structures,  as  divided  into  triangular  elements  for  the’ 

NASTRAN  analyses,  are  presented  in  Fig.  1.  Material  parameters,  were  taken  to  be: 
Youn£f^  modulus , E9  0.1  x 10  dynes/cuT";  density,  p,  2.7  grams/cm'3;  and  Poisson’s 
ratio,  v,  0.31.  (Those  were  experimentally  determined  for  an  aluminum  alloy 
available  in  the  UTRC  shop. ) The  boundary  conditions  for  the  clamped  edges  were 
imposed  by  requiring  that  the  displacements  and  rotations  of  the  grid  points  along 
the  clamped  edges  be  zero. 

Photographs  of  the  experimental  structures,  as  fabricated  to  correspond  to  those 
analyzed  by  NASTRAN,  are  presented  in  Fig.  2.  In  fact,  the  thickness  of  the  experi- 
mental plate  was  0.330  cm  on  the  average. 

The  comparison  of  the  mode  frequencies  of  the  structures,  as  obtained  via 
NASTRAN  and  experimentally,  are  presented  in  Table  I for  the  first  three  vibration 
modes . 

TABLE  I 


STRUClTtRAL  IK) PE  FREQUKNC 1 ES 


Vibration 

Mode 

Plate 

Shell 

NASTRAN 

Experimental 

NASTRAN 

Experimental 

105.6  Hz 

120.6  Hz 

898.7  Hz 

306.0  Hz 

f? 

293.3  Hz 

333.8  Hz 

1280.0  Hz 

1380.0  Hz 

f3 

6 U2.2  Hz 

730.6  Hz 

1*118.0  Hz 

1*080.0  Hz 

The  NASTRAN  computations  of  the  frequencies  for  the  vibration  modes  of  the 
plate  are  all  lower  than  the  experimental  values  by  approximately  lL  percent.  Four 
percent  of  the  difference  is  due  to  the  larger  thickness  of  the  experimental  plate, 
while  the  remaining  ten  percent  is  quite  probably  due  to  a premature  cutoff  of  the 
iterative  solution  for  the  eigenvector  in  the  NASTRAN  computations.  The  consistency 
in  tiie  error  suggests  a repetitive  problem,  and  this  is  supported  by  the  Rayleigh 
quotient  ctiecks  presented  later. 
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The  NASTRAN  oomputat  ion;,  for  tin-  modal  frequenc i os  of 
that  lio  both  above  and  below  the  experimental  value...  It 
first  mode  frequency,  which  in  in  error  by  over  10  percent 
a difficulty  in  adequately  modeling  the  boundary  condition 
attached  to  the  bane. 


t ho  r.hel  1 £ ive  re.  ul  t . 
in  suspect ed  that  the 
, may  be  the  victim  of 
where  the  shell  is 


Holographic  Analysis  and  NASTRAN  Comparison 

Holograms  were*  recorded  of  both  the  plate  fuid  the  shell  while  each  was  vibrat ing 
in  each  of  its  first  three  modes.  In  the  case  of  the  shell,  where  both  transverse 
and  .axial  motions  are  possible,  two  hologram  recordings,  one  with  illumination  from 
the  left  of  the 'shell  and  the  other  with  illumination  from  the  right,  were  made  at 
each  resonance.  The  vibratory  state  was  maintained  exactly  the  same  for  each 
recording  by  means  of  a feedback  circuit  which  employed  an  optical  vibration  probe 
to  detect  the  vibratory  motion,  a bandpass  filter  to  suppress  all  signals  except 
those  in  the  region  of  a resonance,  a limiting  circuit  to  establish  a vibratory 
level,  and  the  conventional  amplifiers  and  electromagnet  driver  to  stimulate  the 
vibration.  When  the  system  was  properly  adjusted,  the  structure's  own  vibration 
mode  would  determine  the  frequency  of  resonance,  tuul  the  limiting  circuit  would 
determine  the  amplitude,  Such  a system  is  inherently  more  stable  than  an  open  loop 
excitation  of  the  resonance  where  a small  drift  of  the  oscillator  frequency  can  pro- 
duct- a significant  change  in  vibration  amplitude.  It  is  important  to  maintain  the 
sane  vibration  level  on  the  structure  for  each  of  the  two  holograms  (with  different 
illuminations)  because  in  their  reconstructions  the  difference  in  fringe  patterns 
between  them  is  used  to  determine  the  motion  transverse  to  the  direction  of  obser- 
vat i on . 


i ". n_t  i • Holographic  Kosult  s 

Reconstruct  ions,  from  holograms  of  the  plate,  recorded  while  it  was  vibrating  h\\ 
each  of  its  first  three  resonance  modes,  are  presented  in  Fig.  s . They  are,  respec- 
tively, the  first  bending,  first  torsion,  and  second  bonding  modes..  From  previous, 
work,  these  patterns  had  been  divided  into  36  sect  ions (a  o x 0 array)  and  the  data 
digitized  and  fitted  to  biquadratic  functions.  For  comparison  to  the  NASTRAN  solu- 
tions., the  latter  were  evaluated  at  the  geometrical  points  corresponding  to  the 
centers  of  the  digitized  sections.  The  differences,  defining  the  mode  fit  between 
the  NASTRAN  and  the  normalized  holographic  results,  are  presented  in  Table  II. 

Pi.  agreement  between  the  two  can  be  seen  to  be  less  than  a few  percent  of  the 
maximum  value  (unity!*  in  all  cases. 

olielj  Holographic  Hy-siil  t s- 


Kecons  t ruct  ions-  from  the  holograms  recorded  of  the  shell  while  it  was  vibrating 
m each  of  its  first  three  modes  are  presented  in  Fig.  U.  These  correspond,  morpho- 
logically, to  the  same  three  modes,  of  the  uniform  plate.  The  photographs  of  these 


l!) 


* 1 ■— 
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Figure  4a  Photographs  of  the  Holographic  Reconstructions  of  the  First  Vibration 
Mode  (806.O  Hz)  of  a 0.317  cm.  Thick  45°  Arc  of  a Cylindrical  Shell 
with  Object  Ream  Illumination  from  the  Left  and  Right  (Clamped 
Edge  at  Bottom) 


io 


* 


reconstructions  were  digitized  using  sections  wlioso  centers  corresponded  to  sixteen 
of  the  NASTRAN  grid  points,  'file  holographic  data  was  processed  to  obtain  vectorial 
displacements  (Ref.  6)  (under  the  assumption  that  motion  in  the  axis  of  the  cylinder 
was  negllbible).  Table  III  shows  the  comparison  between  the  NASTRAN  and  holographic 
results  for  the  first  three  modes  of  the  shell. 


Di f ferences 
in 


TABLE  II 

MODE  FIT  DIFFERENCES:  NASTRAN  VS.  HOLOGRAPHY 


.0005 

-.001*5 

.021*6 

.003b 

-.0010 

. 0202 

.0021* 
-.0019 
. 0206 

.0032 

-.0055 

.0233 

.0027 
. 0061 
. 0224 

. 0031 
.0109 
.0323 

.ool*  1* 

. 0061 

.001*0 

. 0030 

.0023 

-.0011 

.0080 

.000  3 

-.0050 

.0011* 

.0155 

.0150 

. 0200 

.0210 

.0180 

.0190 

.0260 

.0300 

.0060 

. 0060 

.0060 

.0030 

. 0020 

. 00.  ’0 

.001*0 

-.001*0 

-.0101* 

.0161 

. 02  30 

. 0200 

-.0030 

-.0070 

-.0030 

.0010 

.0120 

.0240 

.001*0 

.0030 

.001*0 

. 0080 

. 0000 

. 00  30 

.0010 

-.0120 

-.0296 

.0200 

. 0280 

.0270 

-.0280 

-.0310 

-.0150 

-.021*0 

-.0150 

. 0000 

.0010 

,00.1.0 

. 0020 

.0030 

.0120 

.0010 

.0000 

- . 0190 

-.0110 

.0190 

.0270 

.0140 

-.01*10 

-.0389 

- . 0290 

-.0395 

-.0282 

-.0320 

- . 0020 

.0080 

.0060 

.0010 

- . 00  10 

- . 0060 

.0050 

-.0130 

-.0070 

.0190 

. 0280 

.0190 

.0320 

.0510 

.01*90 

.0270 

.0320 

. 0 320 

System  Checks 

The  next  computations  involved  the  orthogonality  checks  and  Rayleigh  quotient 
checks.  The  former  measure  the  degree  to  which  the  stiffness  functions  between  mode 
are  orthogonal,  and  the  latter  measure  the  ratio  of  maximum  potential  to  maximum 
kinetic  energy  during  the  vibration  cycle  of  a mode  (Ref.  6).  The  orthogonal ity 
parameters  should  be  zero  and  the  Rayleigh  quotients  should  be  unity.  Those 
results  are  presented  in  Table  IV. 
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895 

178 


-.01619 

-.00388 


822 

825 


-.01 303 
-.02961 


.00817 

-.011(36 


-.07816 

.01197 


-.010  30 
-.00388 


-.071+16 

.00751 


-.02628 

.00896 


-.060  13 
-.01111 


-.0071+8  ! 

- . 00  366 

-.01805 

I - . 0 30 12 

.02887 

.00723 

-.01140 

-.02722 

-.02660 

-.02838  j 

-.0 l8s6 

-.02402 

.02735 

- .07566 

-.07800 

.03521 

.07071+ 

. 0 1 388 

- . 00668 

-.0745bl 

.001+35 

-.00302 

-.00279 

.00486 

-.01671+ 

-.02101 

-.0  3 397 

-.0  3050 

.001(71 

.00531 

.00234 

-.009 15 

-.0491+8 

-.0  3 304 

-.04546 

-.05477 

. 12680 

- . 0 49  5 7 

-.04852 

. 1 1050 

.06948 

. 00690 

-.00592 

-.00914 

.00738 

.00514 

. 00022 

.01020 

-.00242 

-.02374 

-.02121 

> -.0  35  82 

-.00114 

1 .00162 

.01154 

.01776 

-.06128 

-.06014 

-.07552 

-.05538 

.0  1647 

-.05390 

-.05400 

. 12700 

-.17110 

.04581 

.01  348 

-.04418 

TABLE  IV 


ORTHOGONALITY  AND  RAYLEIGH  QUOTIENT  CHECK!' 


Mode  Index 
_ (nk) 

Orthogonality  Parameter 

Rayleigh  Quotient 

Plate 

Shell 

Plate 

Shell 

1,  2 

-.01141 

-.000004 

1,  3 

-.00019 

+.085832 

2,  3 

+ .00.109 

.000000 

1,  1 

1.0735 

.0282 

o o 

1.1047 

1.0245 

3,  3 

1.1082 

1.0852 

Structural  Redesign  and  Correlation 

The  next  two  tasks  were  to  redesign  the  structures  and  evaluate  the  results 
via  a second  NASTRAN  analysis.  Both  the  cantilever  plate  and  the  cylindrical 
shell  segment  were  redesigned  subject  to  four  constraints:  a)  lowering  of  the  first 
mode  frequency  by  5 percent  (Awi/ui^  = -.049);  (b)  raising  of  the  second  mode  frequent 
by  5 percent  (A = +.051);  (c)  an  admixture  coefficient  of  the  first  to  the 
second  mode  of  magnitude  0.125  = - .125);  and  (d)  an  admixture  coefficient 

of  the  third  to  the  second  mode  of  magnitude  0.025  (C.,-,  = -.02 !>).  (For  the  plate 

C-  Jj 

the  C.-,^  parameter  was  positive,  and  for  the  shell  it  was  negative.  These  sign 
conventions  were  necessary  to  generate  the  came  shift  in  the  node  line  of  the  second 
mode  on  both  the  plate  and  the  shell  and  they  resulted  only  from  a change  in  the 
sign  conventions  between  ttie  NASTRAN  analyses  of  the  plate  and  the  shell.)  The 
results  provided  structures  which  could  be  fabricated,  and  which  had  root-mean-squaii 
thickness  changes  of  12.7  percent  for  the  plate  and  14.5  percent  for  the  shell. 

The  new  designs,  as  arrays  of  the  fractional  changes  in  thickness  of  the 
triangular  elements  of  the  plate  and  shell,  respectively,  are  presented  in  Figs.  ‘>h 
and  5b . 

Tiie  designs  of  Figs.  5a  and  5b  were  analyzed  by  NASTRAN  to  determine  the  new 
frequencies  and  mode  shapes  of  these  structures.  Tiie  modes  of  the  original  struc- 
tures were  then  used  to  form  a series  to  approximate  (to  least-square-error)  the 
shape  of  tiie  new  second  mode,  and  therefore  to  give  the  resulting  admixture  coeffi- 
cients and  Coj. 

Tiie  final  results  are  presented  in  Table  V,  together  with  the  results  predicted 
by  tiie  more  accurate  perturbation  calculation  mentioned  earlier  in  the  description 
of  the  CHANGNAST  program,  and  considered  further  in  Appendix, 
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Figure  5i>  New  Desigi  for  the  1*5°  Arc  of  a Cylindrical  Shell  as  a Ratio  of 
the  Change  in  Thickness  (Ah/h)  for  Each  of  the  Elements  (Clamped 
Edge  at  Top) 
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TABLK  V 


FINAL  RKSUI/K 


Constraint 


Design 


Plate  Results 


NASTRAN  I Perturb. 


Shell  Results 


NASTRAN  I Perturb, 


Change  in  1st  Mode 
Frequency,  Aio^/u^  -.0it9 


Change  in  2nd  Mode 
Frequency,  Au^/ui,-,  + .0H 


Admixture  Coefficient  I Plate 
1st  Mode  to  2nd  Mode 

C21 

Shell 

-.125 


Admixture  Coefficient 
3rd  Mode  to  2nd  Mode 


-.025 


-.027 


+ .076 


+ .102 


+ . 067 


+ . 03i* 


-.027  I -.028 


The  mode  shape  changes  were  quite  successfully  rendered  by  the  new  design. 
However,  the  frequency  changes  have  fallen  short  of  what  was  required.  All  the 
changes,  however,  arc  in  the  direction  required,  and  for  this  reason  we  my  expect 
that  the  process  would  converge  upon  a desired  design  if  it  were  applied  Herat  vei.v 


SECTION  VI 


fONCLUi'IONf  AND  RKCOMMEN PAT  I ON:' 


Under  the  present  contract,  it  was  sliown  that  it  is  feasible  to  apply  the 
inverse  perturbation  design  method,  originally  developed  at  United  Technologies 
Research  Center  (Refs.  U and  5),  directly  to  the  output  of  a finite  element  com- 
puter analysis  of  a structure.  Thus,  this  new  structural  design  tool,  at  first 
limited  to  a very  simple  geometry  whose  response  characteristics  could  be  suffi- 
ciently characterized  without  the  help  of  a NASTRAN  type  computation,  has  the 
potential  for  application  to  the  mathematically  complex  structure  of,  for  example, 
a turbine  blade.  (The  twist  and  camber  of  such  structures  create  analytical 
difficulties  requiring  the  use  of  finite  element  analysis  to  calculate  their  vibra- 
tion modes  to  the  necessary  accuracy.)  As  described  in  the  previous  section,  the 
initial  step  in  this  direction  has  been  taken  with  the  redesign  of  a shell  shaped 
cantilever  beam.  In  addition,  the  methodology  was  also  applied  to  the  redesign 
of  a flat  cantilever  plate  for  comparison  with  previous  work,  described  in  the 
Appendix  to  the  present  report,  which  only  required  vibration  data  obtained  holo- 
graphically . 

There  is,  of  course,  additional  development,  required  and,  consequently,  it  is 
recommended  that  the  work  be  extended  to  permit  refinement  of  the  analytical 
approach,  together  with  experimental  investigations  concerned  with  the  fabrication 
and  test  of  a redesigned  structure.  Further,  attention  needs  to  be  devoted  to  the 
inclusion  of  additional,  permissable  phy s i cad  design  changes  other  than  the  single 
parameter  of  thickness,  t.o  which  the  technique  has  been  restricted  to  date. 

In  summary,  continued  development  of  the  inverse  perturbation  design  procedure 
would  set  the  stage  for:  a)  redesigning  actual  turbine  blades;  b)  carrying  out 
iterative  design  procedures  without  the  need  for  fabricating  intermediate  proto- 
types; c)  introducing  holographic  data  into  NASTRAN  analyses  to  improve  computer 
modeling  of  turbine  blades  and  other  structures;  and  d)  using  holographic  data 
itself,  rather  than  NASTRAN  output,  for  the  inverse  perturbation  design  of  structure 


APPENDIX 


REDESIGNING  STRUCTURAL  VIBRATION  MODES  BY  INVERSE 
PERTURBATION,  SUBJECT  TO  MINIMAL  CHANGE  THEORY* 


Abstract 

If  first-order  perturbation  theory  for  the  vibration  modes  of  a structure  is  to 
be  inverted  so  as  to  design  a new  structure  having  desired  changes  in  vibration 
modes,  it  is  essential  to  keep  the  structural  changes  as  small  as  possible.  This 
not  only  helps  keep  the  changes  within  the  approximation  of  first-order  perturbation 
analysis,  but  also  avoids  ill-conditioning  of  the  simultaneous  equations  that  must 
be  solved  for  a new  design.  Calculus  of  variations  may  be  used  to  find  the  minimal 
changes,  and  specifies  a unique  set  of  functions,  derived  from  pairs  of  vibration 
modes,  that  are  optimal  for  expressing  the  structural  changes  needed  to  accomplish 
a given  set  of  changes  in  vibration  modes.  This  theory  is  put  into  practice  for  a 
cantilever  plate,  for  which  vibration-mode  data  is  obtained  by  holographic  vibration 
analysis . 

Introduction 

Two  previous  papers  (Refs.  U and  5)  have  outlined  both  a general  first-order 
perturbation  theory  for  the  vibration  modes  of  a structure  and  the  inversion  of  that 
theory  for  the  purpose  of  structural  redesign.  Specifically,  the  first  answered 
the  question:  How  will  a small  change  to  a structure,  whose  vibration  modes  are 
known,  affect  those  vibration  modes?  The  second  answered  the  question:  How  should 
a structure,  whose  vibration  modes  are  known,  be  modified  to  alter  certain  modes? 

The  use  of  first-order  perturbation  theory  to  answer  the  first  question  may  not 
be  practical  because  finite-element  computer  analyses  can  often  provide  such  answers 
more  accurately.  This  is  because  first-order  approximations  often  become  more 
error  prone  in  vibration  theory  than  modeling  approximations  by  finite  elements. 

With  the  second  question,  however,  perturbation  methods  may  have  an  advantage  if 
a considerable  number  of  simultaneous  changes  to  vibration  modes  are  required. 

This  is  because  the  computation  time  required  for  finite-element  analyses  can  make 
hunt-and-try  methods  impractical.  Still,  dependence  upon  first-order  variational 
theory  of  vibration  modes  can  make  an  inverse-perturbation  design  procedure  quite 
inaccurate  unless  the  structural  changes  are  kept  small.  A method  is  needed, 
therefore,  to  find  the  smallest  physical  change  to  a structure  that  will  alter 
its  vibration  modes  in  the  desired  manner. 


* This  appendix  was  authored  by  K.  A.  Stetson,  I.  R.  Harrison  and  G.  E.  Palma. 
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The  next  section  of  this  appendix  presents  an  application  of  the  calculus  of 
variations  to  the  problem  of  finding  the  minimum  structural  change  necessary  to 
achieve  a desired  set  of  changes  in  the  vibration  modes.  The  remainder  of  the  appen- 
dix deals  with  the  application  of  this  theory  to  the  redesign  of  a cantilever 
plate,  based  upon  data  acquired  by  holographic  vibration  analysis.  The  design 
parameter  permitted  to  vary  was  the  thickness  of  the  plate,  although,  in  principle, 
the  length  and  breadth  of  the  plate  could  be  used  as  design  variables  also.  This 
application  required  the  development  of  both:  _l)  the  necessary  computer  program 
to  perform  the  inverse  perturbation  design  routine;  and  2)  a program  to  analyze 
holographic  interferograms  of  vibration  modes  so  as  to  obtain  accurate  second 
derivatives  of  the  normal  mode  functions.  These  latter  were  required  to  model  the 
changes  in  bending  strain  due  to  changes  in  plate  thickness.  In  the  end,  new 
designs  were  actually  generated  and  built  as  cantilever  plates  of  nonuniform  thick- 
ness . These  new  plates  were  then  analyzed  to  compare  the  actual  changes  in  normal 
vibration  modes  with  the  changes  prescribed. 


Theory 


Minimization  of  Inverse  Perturbation 


The  perturbation  theory,  as  presented  in  Ref.  4 can  be  summarized  as  follows. 
Let  us  assume  that  a small,  distributed  change,  f (x) , is  made  to  a structure. 

For  every  pair  of  vibration  modes  of  the  structure,  $n(x)  and  $^(x),  there  exists 
a function,  0nk(x)>  which  describes  the  effect  the  change  will  have  on  these  modes 
by  the  integral 

a * fLenk(*)  f(x)  dx- 

°nk  J nk  (1A) 

0 

The  variable  x is  taken  to  represent  one,  two  or  three  dimensions  depending  upon 
the  nature  of  the  structure,  and  L and  0 denote  its  boundaries.  The  parameter, 
Ank,  describes  change  in  either  the  square  of  the  natural  frequency  if  n=k,  or 
admixture  between  modes  if  n^k,  i.e., 


snn  " V = cnk("n-“k)/“n“k,  <“> 

where  u and  w.  are  the  natural  frequencies  (in  radians)  of  the  nth  and  kth  modes, 
n k 

and  C , is  the  admixture  coefficient  that  describes  the  amount  of  the  modeshape 
nk  k 

that  the  4>n  mode  subsumes  as  a result  of  the  structural  change.  The  mode  admixutre 

is  described  by 
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«;(x)  . ♦„(*>  * £ Cnk»k(x) 


( 3A ) 


where  £ (x)  is  the  new  shape  of  mode  I , ami  k=n  is  omitted  from  the  summation,  The 
. . n , . n 

Junctions,  Oj^tx;,  which  may  be  called  perturbation  functions,  need  not  be  explicitly 
defined  at  this  point. 


In  redesigning  the  structure,  we  wish  to  specify  certain  changes  in  modal 
frequencies  Au  , and  certain  changes  in  mode  shapes  expressable  by  admixture  cc 


cients  C . ; these  changes  specify'  certain  A's  via 


(2A).  The  goal,  of  course. 


is  to  find  a suitable  function  f(x),  which  will  yield  these  A's  via  Eq.  (1A)*. 
Unfortunately,  there  may  exist  a wide  clioioe  of  functions  that  all  yield  these  A's. 
Let  us  seek,  however,  to  find  that  function  which  will  make  the  minimum  change  to 
the  structure.  If  we  follow  the  procedure  outlined  by  Schultz  and  Melsa  (Ref.  7), 
we  first  require  a suitable  index  of  the  change  to  the  structure;  for  example, 
the  integral  of  the  square  of  f(x).  This  may  be  referred  to  as  a functional  of  the 
function,  f(x),  in  that  for  any  given  function,  the  functional  takes  on  a given 
value.  If  the  variation  of  the  functional  exists  and  the  functional  takes  on  a 
minimum  for  a given  function,  then  the  variation  of  the  functional  will  be  zero  for 
that  function.  Wo  seek,  therefore,  a function  that  will  satisfy  the  condition 


6 JV(x)  dx 
0 


However,  we  also  want  the  function  to  yield  a given  set  of  A’s  via  a given  s 
Kqs . (1A),  i.e.,  when  integrated  times  a given  set  of  perturbation  functions 
Onk(x).  These  conditions  are  contraints  upon  our  choice  of  functions,  f(x). 
may  be  introduced  into  Eq.  (1*A)  by  the  technique  of  undetermined  multipliers 
us  use  the  index  p for  the  node  pair  nk,  and  express  the  constraints-  i:.  the 
form: 


They 

Let 


fol  lev.  : ng 


0 = 


©p(x)f(x) ]dx 


(SO 


Let  us  assume  that  we  have  specified  a set 
of  tiie  structure.  Since  eacli  of  the  P Eqs 
by  an  arbitrary  factor,  2X^ , and  add  it  to 
of  the  functional.  This  gives 


of  E constraints  as  our  goal  for  redesign 
(5A)  equals  zero , we  may  multiply  each 
Eq,  (1.A)  without  changing  the  variance 


6 (f2(x)  + 2lXp[Ap/L 


0p( x) f ( x) ] } dx  = 0. 


l oA ' 


r 


f 


T 


1 


Since  the  variational  equation  is  in  the  form 

6 /L  V (f(x)  ,X)  dx  = 0, 
0 N 

we  are  led  to  the  Euler  equation 


3 y/3f  » 0. 


(7A) 


This  leads  to  the  solution  for  the  desired  function,  f(x). 


11 

1 


f(x) 


Vx)- 


(8a) 


The  perturbation  functions  take  on  an  interesting  role  in  this  problem  in  that  they 
not  only  may  be  used  to  determine  the  effects  of  any  structural  changes,  they  them- 
selves form  a set  of  optimal  functions  for  formulating  changes  to  the  structure. 

To  determine  the  coefficients  of  the  series  in  Eq.  (8A),  the  A^'s,  substitute 
for  f(x)  in  Eq.  (1A).  ^ 


A 


P 


P 

l \n 

1 P 


/ 


L 

0pU)  0q(x) 


dx. 


(9A) 


There  are  P such  equations  each  of  which  has  F unknown  X's,  for  which  we  may  solve 
by  the  matrix  equation 


(x)  = [*]  1 (A), 


(10A) 


where  the  elements  of  H’ , which  may  be  called  the  perturbation  matrix,  sire 


/L0n(x)  0o(x)  dx. 
0 M M 


(HA) 


(Parentheses  denote  vectors 
The  final  structiu'al  change 


or  single  column  matrices  and  square  brackets  matrices.) 
may  be  expressed  by 


i 


1 


A 
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f ( X ) = (a)T[H  ]-  1 (o(x)).  ) 

II’  the  perturbation  functions  are  independent  over  the  boundaries  of  the  structure, 
it  can  be  shown  (Ref.  3)  that  the  matrix  generated  by  Eq.  (llA)  must  be  nonsingular. 
In  general,  tiie  perturbation  functions  are  independent  even  though  they  are  not 
orthogonal.  This  is  true  even  if  complementary  pairs  of  perturbation  functions , 

0nk  and  0^,  , are  included,  because  these  functions  are  not  proportional  unless  run 
mass  or  pure  stiffness  changes  are  made.  The  rank  of  the  perturbation  matrix  is 
determined  solely  by  the  number  of  constraints,  or  redesign  changes,  required,  and 
not  by  tiie  number  of  degrees  of  freedom  of  the  structure.  This,  generally  means 
that  the  perturbation  matrix  will  not  be  very  large. 

Formulating  the  problem,  as  above,  allows  one  further  insight  into  tiie  process 
of  structural  redesign.  It  is  quite  easy  to  show  that  the  absolute  minimum  root  - 
mean-square  change  to  tiie  structure  needed  to  accomplish  a suggested  redesign  is 

Min  {.f  f2(x)  dx}1/2  - {(a)T[*]' H a)  f /2 

This  gives-  a measure  of  the  difficulty  involved  in  making  any  suggested  change  in 
the  vibration  modes  of  a structure.  It  will  not,  however,  indicate  any  s.pecific 
violations  of  physical  reality,  such  as  negative  plate  thicknesses. 

Application  to  Plate  Design 

For  an  isotropic  plate,  whoso  thickness  we  may  vary,  we  nay  (as  in  '.‘of.  •») 
write  tiie  following  express  ion  for  the  change  parameter,  Ank, 

Ank  = (v,kMk)"1JJC(*k)T^  (*n)((1+Ah/h)3"1)_"nph*nVh/h]  dxdy>  UJ,A' 


M.  i.-  the  modal  mass  of  the  kth  mode,  defined  by 
K.  1 


Mk  = / ptuk  dxdy  ’ 


(i  s\) 


where  p is  the  volume  density  of  the  material,  and  h is  the  slate  thickness. 
Eq . (lUA),  All  is  the  change  in  plate  thickness  and  (^n)  is.  defined  by 


. n 


1 1 


(16A) 


where  superscripts  denote  partial  differentiation  with  respect  to  the  superscript 
variables,  and  G is  defined  by 


G 


1 v 0 
v 1 0 

0 0 2 (1-v) 


Eh3/12(l-v2)  , 


(17A) 


„ T 

where  K is  Young's  modulus  and  v is  Poisson's  ration.  (4>^)  denotes  the  transpose 
of  the  column  matrix,  formed  for  the  kth  mode,  as  defined  in  Eq.  (l6A).  The  per- 
centage change  in  frequency  of  a mode  and  the  admixture  coefficients  between  modes 
may  be  calculated  from  the  change  parameters  by 


A“n/a,n  = (1+  ann) 1/2-1  ‘ and  (i8a) 

Cnk  = Ank“n“k/K-“k)  ‘ (iqa) 


I, 


‘ 1 


I 


To  put  Eq.  (lUA)  into  the  form  of  Eq.  (lA),  we  must  obtain  Ah/h  (the  percentage 
change  in  thickness)  as  a common  factor  of  both  terms  in  the  integrand.  This  can 
be  done  by  taking  the  first-order  approximation 

((l+Ah/h)3-!)  = 3Ah/h.  (20A) 


From  this,  the  perturbation  function  becomes 

°nk  = K“kV_1  [3(*n)T-G  (*k)-u)n  phVk]> 


This  is  tlie  function  that  may  be  used  in  Eqs.  (11A)  and  (12A)  with  Ah/li  replacing 
f(x). 
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The  use  of  approximation  (20A)  yields  a design  procedure  that  makes  design 
changes  that  accurately  model  changes  in  mass,  but  only  approximately  model  change:-, 
in  stiffness.  Alternatively,  we  could  choose  linear  changes  in  the  cube  of  thick- 
ness, A(h^)/h\  as  our  structural  change , f(x),  and  perform  the  approximation  in 
the  second  term  of  the  integrand  of  Eq.  (l^A).  This  leads  to  a perturbation  function 
that  is  simply  1/3  the  value  of  that  in  Eq.  (21A).  Thus,  we  may  consider  either  of 
two  possibilities  for  the  redesign  of  an  isotropic  plate  based  on  exact  changes  in 

mac.;'.,  (Ah/h)  , or  on  exact  changes  in  stiffness,  ( Ah/h ) The  two  are  related  by 
m s 

( Ah/h) s = (1  + 3(Ah/h)J1/3  . (22A) 

Ideally,  it  would  be  desirable  to  make  design  changes  that  were  modeled  exactly  in 
both  mass  :uid  stiffness.  Attempting  to  do  this  leads  to  an  Euler  equation  that 
cannot.be  solved  meaningfully  in  the  simple  form  of  Eq.  (8A).  Nonetheless,  once  a 
design  has  been  obtained,  Eq.  (lUA)  can  be  used  to  generate  a more  exact  estimate 
of  the  values  of  Ank  to  be  expected. 


Experimental  Application 


wenernl  Cons i derations 


In  choosing  an  experimental  application  of  the  proposed  design  technique,  we 
considered  the  following:  There  was  considerable  interest  in  a method  that  could 
function  with  experimentally  determined  data,  specifically,  holographically  obtained 
data.  There  was  interest  in  structures  that  resembled  turbine  blades;  however,  for 
simplicity,  we  preferred  to  work  with  a plate  structure.  This  led  to  the  choice  of 
a cantilever  plate.  Holographic  vibration  analysis  would  provide  displacement  data 
for  the  vibration  modes  which  would  require  numerical  processing  to  obtain  second 
derivatives.  Errors  in  displacement  data  generally  lead  to  severe  errors,  in  compu- 
tation of  second  partial  derivatives.  Since  our  second  partial  derivative  functions 
were  to  be  integrated,  to  form  a perturbation  matrix,  or  to  model  distributed 
changes  in  the  thickness  of  the  plate,  the  average  value  of  the  second  partial:',  over 
segments  of  the  plate  would  be  more  important  than  their  peak  values.  (The  latter 
are  more  important  for  prediction  of  failure  due  to  stress.)  This  led  to  the  idea 
of  dividing  the  plate  into  sections,  fitting  the  data  within  each  section  to  a bi- 
quadratic function  by  least-squure-error  theory,  and  differentiating  this  function 
to  obtain  the  average  value  of  the  second  derivative.  Furthermore , there  would  be 
similarities  between  this  approach  and  the  approach  we  would  have  to  take  if  wo 
were  to  use  vibration  modes  computed  by  finite  element  analyses  as  our  input. 


MODES  Program 


A computer  program,  named  MODES,  was  written  to  obtain  the  second  partial 
derivatives  of  the  modes  of  a plate.  The  program  accepts  data,  in  the  form  of  x 
and  y coordinates  and  a fringe  number,  from  photographs  of  holographic  reconstruction 
of  the  mode  functions.  It  then  fits  this  mode  data  for  least-square  error  to  a 
quadratic  function  of  two  coordinates. 

Because  the  second  partial  derivatives  are  required  for  the  computation  of  the 
perturbation  functions,  the  mode  shape  must  fit,  to  lowest  order,  a quadrat i c func- 
tion in  two  dimensions.  A bi-quadratic,  rather  than  a bi-cubie  fit,  is  used  in 
order  to  simplify  the  subsequent  calculations,  and  to  minimize  the  data  error.  A 
cubic  fit  would  follow  the  data  errors  more  closely  because  of  its  added  degrees 
of  freedom.  The  bi-quadratic  function  used  may  be  expressed  as 

<j>  = Ax2  + By2  + Cxy  + Dx  + Ey  + F,  (13A) 


where  <1>  is  the  mode  value  at  any  point  x,  y and  2A  is  the  second  partial  derivative 
with  respect  fo  x,  |J>XX;  2B  is  the  second  partial  derivative  with  respect  to  y,  ; 
and  C is  the  cross  partial  derivative,  4>xy.  D and  E are  the  average  first-partial 
derivatives  with  respect  to  x and  y,  respectively,  and  F is  the  average  value  over 
the  section. 

We  are  approximating  a mode  function  of  a structure  which  has  continuous  and 
smooth  first  and  second  partial  derivatives.  Whereas  the  vibration  amplitude 
represented  by  a fringe  is  exact,  there  is  error  encountered  in  locating  the  center 
of  a fringe.  In  order  to  minimize  the  effect  of  this  error,  we  must  average  over 
the  data.  In  this  pirocess , the  fewer  the  degrees  of  freedom,  the  more  the  error 
effect  is  minimized,  which  is  wliy  the  data  is  fit  to  a bi-quadratic,  which  has  con- 
stant second  partial  derivatives.  Because  the  fit  to  the  vibration  data  yields 
constant  second  partial  derivatives  over  the  area  considered,  it  is  necessary 
to  analyze  the  structure  in  sections,  which  gives  an  average  fit  for  each  section 
of  the  mode.  The  constraint  that  the  total  function  be  continuous  between  sections 
could  be  added,  and  a solution  obtained;  however,  doing  this  makes  the  computation 
more  complex  since  the  fits  for  all  sections  must  be  solved  simultaneously.  For 
simplicity,  data  from  each  section  of  the  structure  was  fit  independently,  and 
continuity  sacrificed.  To  get  a good  approximation  to  the  mode  function  without 
the  continuous  and  smooth  boundary  constraints  between  sections,  many  sections 
should  be  used.  We  must  strike  a balance,  however,  between  getting  a good,  smooth 
overall  fit  to  the  mode,  and  getting  enough  data  per  section  to  minimize  data 
error.  This  was  done  empirically  by  recording  holograpms  with  high  fringe  den- 
sities and  dividing  the  structure  so  that  at  least  fifty  data  points  could  be 
taken  per  section. 


After  the  data  from  each  section  of  the  strucutre  has  been  fit  to  a bi-quadrut 
function,  a computer  storage  file  is  formed  for  each  mode,  with  the  coefficients 
of  t!ie  particular  mode  fit  and  the  center  coordinates  for  each  section  of  the  struo 
ture  stored  in  that  file. 

In  actual  use,  the  program  was  run  on  a time-sharing  computer  terminal  ami  a 
mode  file  name  and  fringe-order  type  were  chosen.  For  each  section,  data  was 
entered  as  an  x-y  position  and  fringe  number.  A coordinate  shift  was  made  on  cacti 
set  of  data  points,  relative  to  the  center  of  the  section.  Each  data  point,  and 
its  corresponding  fringe-order  value,  were  used  to  compile  matrix  coefficients, 
the  contributions  to  which  were  linearly  additive  for  each  data  point.  After 
adding  the  term  for  that  point  into  the  matrix,  each  data  point  is  discarded  as 
the  next  point  is.  read.  The  matrix  and  the  vector,  formed  and  solved  for  each 
sect  ion  are 
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where  the  letters  A through  V are  coefficients  of  the  fitted  polynomial  in  Kq.  ( 1 A ) 
xn  and  y are  the  coordinates  of  the  data  point  of  the  fringe  whose  order  generates- 
the  value  Sln , and  N is  the  total  number  of  data  points.  After  sufficient  data 
has  been  entered  for  a section,  the  user  types  END,  which  signals  the  computer  to 
solve  for  the  coefficients  A through  K of  <l>n.  This  process  is  repeated  for  each 
section  of  the  structure,  and  all  sections  are  included  in  the  mode  fit  file. 

Test. ing  the  MOPED  Program 

The  MOPED  program  was  tested  for  its  accuracy  in  calculating  second  derivative 
and  its  sensitivity  to  data  error  by  taking  data  from  a known  mode  function  which 
was  not  of  quadratic  form.  The  data  was.  generated  numerically  from  an  analytical 
function,  and  these  values  were  used  in  the  mode  fitting  program.  Then,  random 
error  was  added  to  the  data  and  the  sensitivity  to  the  error  was  determined. 


The  theoretical  solution  for  a plate,  hinged  ut  all  edges,  war.  used  for  this 
test.  The  two  vibration  modes  were 


$1  = sin(x)  sin(y) , and 
$2  = sin(x)  sin(2y)  . 


(p‘>A) 


Contour  plots  of  these  modes  are  shown  in  Fig.  1A.  The  first  quadrant  of  this 
plate  was  divided  into  25  sections,  mid  the  data  from  these  sections  were  used  in 
the  MODES  program.  The  exact  values  of  the  second  partial  derivatives,  averaged 
over  a section,  and  the  exact  values,  calculated  at  the  center  point,  were  com- 
puted from  the  analytical  equations.  Kundom  error  was  then  added  to  the  same  data 
mid  used  in  the  MODES  program.  (The  error  was  added  to  the  fringe  values.,  'although 
the  actual  error  occurs  in  the  x-y  values,  because  the  real  error  occurs  as  a per- 
centage of  the  fringe  spacing  and  varies  in  sine  from  fringe  to  fringe.  Do  long 
as  the  error  is  small,  this  can  be  approximated  by  an  error  of  uniform  sine  in  the 
fringe  values.)  Different  percents  of  error  were  tested,  and  the  results  (s)iown 
in  Table  IA)  indicate  that  less  than  5 percent  error  is  desirable. 


TABLE  TA 


RANDOM  ERROR  EKKKCTB  ON  MOPED.  BROOKAM 


Sect.  1 

1 percent 

5 percent 

10  percent 

AXX 

<9  l 

-.0845 

-.0877 

-.0321 

*yy 

-.1139 

-.0696 

-.349 

»xy 

+.9077 

+.8955 

+.8589 

Sect.  2 

XX 

♦ 

-.1843 

-.2581 

-.3508 

*yy 

-.1842 

-.1135 

-.0957 

*xy 

+.7247 

+.7261 

+.5735 

Sect.  3 

*xx 

-.0715 

-.0365 

-.1099 

♦yy 

-.0695 

-.0558 

-.1479 

,xy 


+.6188 


+.6058 


+.6636 


The  repeatability  of  the  MODES  program  when  a section  is  redigitized  is  also 
a measure  of  its  sensitivity  to  random  error.  A photograph  of  a holographic 
vibration  mode  was  placed  on  the  digitizer  table,  and  the  same  region  was  digitized 
several  times.  At  most  there  was  a two  percent  difference,  based  on  successive 
digitizations,  in  the  perturbation  function  computed  for  the  segment. 

To  compare  ttie  mode  fit  for  data  from  an  actual  structure,  with  an  analytical 
solution,  we  used  a disk  fixed  in  the  center.  The  analytical  solution  for  a circular 
disk  with  free  edges  is  available.  A photograph  of  the  holographic  reconstruction  of 
a two  diameter  mode  was  digitized  and  fitted  in  each  of  three  sections.  The  sections 
were  at  different  radii  and  angles  from  the  node  line  (see  Fig.  2A).  The  calculated 
perturbation  function  for  each  section  was  compared  with  the  exact  analytical 
value  at  the  center  point  of  each  section.  The  percent  differences  for  each 
section  were:  section  1,  12.5  percent;  section  2,  10. U percent;  and  section  3, 

26.0  percent. 

Each  calculated  value  is  for  the  point  witli  the  specified  radius  and  angle,  and 
ttie  mode  fit  is  an  average  over  an  area.  If  the  area  is  decreased  in  size,  the  fit 
should  approach  ttie  exact  solution. 

Cl  LAMP ED  Program 

Another  program,  CHANGES,  solves  the  equations  described  in  the  section  on 
Theory  to  calculate  the  thickness  changes  required  in  a structure  to  make  desired 
changes  in  the  vibration  modes.  The  user  specifies  which  modes  are  to  be  changed, 
and  by  what  amount.  The  change  parameters  specified  (i.e.,  A-q,  A-^-,,  etc.)  deter- 
mine which  perturbation  functions  are  required.  The  program  accepts  the  mode 
coefficients  (and  normalizes  them),  the  plate  thicknesses,  and  other  constants 
needed  to  form  each  perturbation  function.  For  simplicity,  each  perturbation  func- 
tion has  a constant  value  for  each  section.  This  requires  that  the  modal  component 
of  the  perturbation  function,  $ 4>k,  be  averaged  over  each  section.  As  a result, 
each  section  will  be  assigned  a constant  change  in  thickness,  and  the  new  plate 
will  have  a thickness  profile  with  steps  of  differing  height.  Using  the  perturba- 
tion functions,  the  computer  forms  the  perturbation  matrix,  and  from  the  change 
parameters,  the  vector  of  A's.  It  then  solves  for  the  new  thicknesses  of  a plate 
having  the  desired  mode  changes,  based  upon  exact  changes  in  mass,  (Ah/h),n,  and 
prints  out  the  RMS  thickness  change.  The  program  then  goes  back,  and,  using  the 
new  plate  thicknesses,  computes  all  the  new  mode  changes;  even  those  where  change 
was  not  specified,  via  the  more  accurate  method  of  Eq.  (lUA). 

The  perturbation  functions  for  the  simply  supported  plate  were  calculated 
analytically  on  a programmable  calculator.  These  perturbation  functions  were  cal- 
culated for  the  two  mode  functions  in  Eqs.  (3A)  and  (1(A),  and  the  perturbation  func- 
tions can  be  seen  in  Fig.  3A  (a,b,c,  and  d).  (Note  that  0-^  and  Cb^,  for  example, 
have  their  maximum  values  at  the  corners  of  the  plate  even  though  the  mode  functions 


38 


(3.14.0) 


(3.14.01 


are  zero  there.  This  results  from  the  high  bending  that  occurs  in  the  corners.) 
next,  the  mode  coefficients  for  one  quadrant  of  the  simply  supported  plate,  as 
calculated  from  computer-generated  data  points  (see  MODES  sections),  were  used  in 
the  pert urbat ion  function  part  of  the  CHANGES  program.  These  perturbation  function 
values  were  compared  with  computer  calculated  values  for  the  exact  analytical  solu- 
tion of  the  perturbation  function  for  the  plate.  The  results  are  shown  in  Fig. 

**A  (a,b,c,  and  d).  Tabulated  are  the  exact  values  at  the  center  of  the  section 
and  the  percentage  error  of  the  numerically  calculated  values . The  results  indicate 
that  a large  percentage  of  error  occurs  only  where  the  perturbation  function  is 
small. 


In  order  to  check  if  changes  in  thickness  gave  the  correct  mode  perturbations, 
the  thicknesses  of  two  corners  of  the  simply  supported  plate  were  thinned  by  30  per- 
cent. Again,  the  mode  coef ficienets  for  the  simply  supported  plate,  as  calculated 
from  computer-generated  points,  were  used  in  the  perturbation  function.  This  time 
the  mode  perturbations  were  compared  with  the  mode  perturbations  calculated  analyt- 
ically. The  results  are  in  Table  II  A. 


TABLE  I IA 

PERTURBATION  COMPARISONS 


— 

Computed 

Analytic 

Percent  Diff. 

“1 

151.779  Hz 

151.645  Hz 

+.08  pet. 

w2 

387.233  Hz 

386.378  Hz 

+.22  pet. 

C12 

+.012923 

+.013311 

-.291  pet. 

C21 

-.01166 

-.01277 

-8.69  pet. 

It  was  concluded  that  the  CHANGES  program  worked  properly,  and  that  it  and  the 
MODES  program  should  be  applied  to  an  experimental  test. 

Implement  at ion 

The  structure  chosen  for  redesign  was  a 12.7  cm  by  15 .2U  cm  cantilever  aluminum 
plate,  clamped  on  one  end,  and  excited  by  electromagnetic  transducers  stimulating 
eddy  currents  in  the  aluminum.  The  excitation  frequency  was  adjusted  until  the 
plate  was  vibrating  in  one  of  its  normal  modes,  identified  by  observation  through  a 
speckle  interferometer.  A suitable  amplitude  was  set,  and  holograms  were  recorded 
for  the  first  five  modes,  the  first  three  of  which  were  presented  as  Fig.  3 in  the 
main  body  of  the  report.  The  frequency  of  the  maximum  response  for  eacli  mode  was 
recorded  (fx  = 120.6  Hz,  f0  = 333.3  iiz,  f,  = 730.6  Hz)  and  photographs  were  taken 
of  each  holographic  reconstruction  and  enlarged  prints  made. 


Fringe  data  from  the  mode  photographs  was  extracted  by  means  of  an  x-y 
position  digitizing  table  (279  ram  square)  made  by  Summagraphics , Inc.,  which  was 
accurate  to  0.005  inches  (0.13  mm).  This  consisted  of  an  active  table,  a position 
display,  and  a cursor,  which  itself  consisted  of  cross  hairs,  position  button,  and 
flag  buttons.  The  cross  hairs  were  centered  over  a point  on  the  table  and  the 
position  button  pushed.  The  X and  Y values  (in  left-hand  coordinates)  for  that 
point  on  the  table,  then  appeared  on  the  display,  and  were  sent  to  the  computer. 

An  overlay  grid  was  drawn  which  had  the  same  perimeter  as  the  plate  in  the 
photographs  and  was  divided  into  36  sections  (6  x 6),  with  the  center  of  each 
section  marked.  This  overlay  grid,  positioned  and  permanently  fixed  to  the 
digitizer  table,  was  used  for  all  three  modes  to  delineate  the  sections  of  evalua- 
tion. Each  mode  photograph  was  slid,  in  turn,  under  the  grid  for  evaluation, 
insuring  that  the  orientations,  sizes,  positions,  and  centers  of  each  section  would 
be  tiie  same  for  each  mode. 

For  each  section,  tiie  center  point  was  recorded  and  the  fringes  were  digitized, 
one  at  a time.  An  initial  fringe  number  (tiie  lowest)  was  entered  into  tiie  computer, 
and  after  digitizing  eacli  fringe,  the  fringe  number  could  be  increased  by  one  for 
tiie  next  fringe  by  pushing  the  cursor's  flag  button,  or  by  entering  it  separaetly 
via  the  terminal.  After  eacti  section  was  finished,  an  instruction  (END)  initiated 
computation  of  tiie  mode  coefficients  for  that  section,  and  the  process  was.  repeated 
for  the  next  section. 

After  the  mode  coefficients  for  all  sections  of  all  three  modes  were  found, 
the  A,  B,  and  C coefficients  for  the  first  three  modes  were  plotted  on  isometric 
paper.  These  drawings  were  examined  for  smoothness,  and  sections  that  appeared 
irregular  were  redigitized  to  ciieck  their  accuracy.  If  new  values  were  more  regular 
they  were  used  instead. 

The  exact  thickness  of  eacti  section  of  the  plate  was  measured  (the  thickness  of 
the  plate  was  roughly  .3175  cm,  but  the  programs  allow  for  different  thicknesses 
for  eacli  section),  and  the  other  constants  determined  experimentally  (E  = 6.1  x 10 
((jynes/cm"  ) , p = 2.7  (g/cm^),  and  v = .31).  The  mode  coefficients,  plate  thickness, 
arid  other  constants  were  used  in  tiie  CHANGES  program.  Several  different  designs, 
with  different  mode  perturbations  were  examined,  and  one  sucli  design  was  chosen 
for  fabrication  via  the  (Ah/ii)  model.  The  design  chosen  had  change  parameters 
ALi  = -.10,  A00  = +.10,  A21  = -.5991*,  and  A„3  = .0380.  Tiie  first  two  correspond 
to  changes  in  frequency  of  the  first  two  modes  of  Aw^/io^  = -.106  and  Aw^/ui.,  = +.095* 
Tiie  second  two  correspond  to  admixture  coefficients  of  = -.21*91  and  C.,_,  = -.050b 
which  displace  the  node  line  of  the  second  mode  as  shown  in  Fig.  5A.  The  ratios 
Aii/ii  for  eacli  section  are  presented  in  Table  TTTA,  arranged  in  the  6 x 6 format  of 
the  sections  on  tiie  pilate.  Tiie  RMS  change  in  thickness  was  0.281*. 
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TABLE  TTTA 


FRACTIONAL  THICKNESS  CHANCES 


Clamped  End 


-.109 

-.169 

-.356 

-.294 

-.084 

+ .196 

+ .011 

+ .028 

-.007 

0 

1 

+ .003 

+ .016 

1 

O 

cn 

00 

+ .111 

+ .193 

+.222 

+ .147 

+ .047 

-.230 

+ .006 

+ .135 

+ .294 

+ .265 

+ .218 

-.506 

-.107 

+ .135 

+ .284 

+ .404 

+ .281 

-.847 

-.352 

+ .153 

1 

+ .506 

+ .504 

Free  End 


Figure  6A  shows  the  new  plate,  cemented  and  bolted  to  the  same  base  block  used 
for  the  uniform  plate,  and  located  in  position  with  respect  to  the  hologram 
recording  apparatus.  The  frequencies  of  the  first  three  modes  for  this  plate  were 
106.3  Hz,  381.5  Hz,  and  770.6  Hz,  respectively,  and  a photograph  of  the  holographic 
reconstruction  of  the  new  second  mode  is  presented  in  Fig.  7-A.  The  admixture 
coefficients  were  determined  by  digitizing  section  3*t , the  free  end  section  that 
contained  the  displaced  node  line,  and  processing  the  data  with  the  MODES  program 
to  get  the  bi-quadratic  fit  coefficients.  These  were  then  normalized  and  used 
with  the  normalized  bi-quadratic  fit  coefficients  for  that  section  from  the  first 
three  modes  of  the  original  uniform  plate  to  form  a set  of  simultaneous  equations 
expressing  the  new  mode  shape  as  a linear  combination  of  the  old  mode  shapes. 
Solution  of  these  equations  gave  the  admixture  coefficients  = -0.217  and 

c23  = -0.0935. 


Given  the  plate  with  modified  thickness,  as  tabulated  in  Table  TTA,  there  are 
three  ways  to  calculate  the  expected  perturbations  of  the  vibration  modes.  First, 
we  may  integrate  Ah/h  times  0 , as  defined  in  Eq.  (21A).  Second,  we  may  compute 
the  corresponding  variable  A(iP)/h^  = ( (l+Ah/h ) 'l-l ) and  integrate  its  product  with 
®nk/3*  (The  latter  will  correctly  model  changes  in  stiffness,  whereas  the 
former  will  correctly  model  changes  in  mass.)  Third,  we  may  use  both  parameters 
in  Eq.  (lUA)  to  get  the  most  accurate  estimate  of  the  first-order  mode  perturbations 
Tiie  results  of  these  three  methods  are  tabulated  in  Table  TVA,  together  with  the 
experimentally  determined  values.  Also  tabulated  are  the  root -mean-square  differ- 
ences between  the  experimental  and  eacli  of  the  calculated  values. 


51 


(Note  that  the  perturbation  of  the  frequency  of  the  third  mode  has  been  included, 
even  though  it  was  a parameter  whose  change  was  not  specified.  The  perturbations 
calculated  via  Ah/h  (except  for  Aw^/uO  were  the  design  objectives,  and  the  RMS 
difference  between  these  and  the  experimentally  obtained  results  was  .0365.  Tins 
is  15  percent  of  the  root-mean-square  of  the  design  changes.  The  last  column 
presents  the  absolute  difference  between  the  originul  design  requirements  and  the 
experimentally  obtained  values.) 


Considering  the  complexity  of  the  design  requirements,  and  the  size  of  the 
changes  specified  (admixture  between  modes  1 and  1 was  25  percent),  we  should 
regard  the  results  as  quite  favorable.  All  constrained  parameters  did  change  in 
the  correct  direction,  and  it  is  reasonable  to  assume  that  a second  iteration  of 
the  method,  based  on  data  taken  from  the  second  plate,  would  result  in  a plate  whose 
vibratory'  performance  would  be  quite  close  to  specification.  It  should  be  noted 
that  ttie  largest  change  parameter,  C01  , seems  to  have  dominated  the  design  at  this 
step.  The  last  column  of  Table  TVA,  which  would  provide  the  input  requirements 
for  the  next  iteration,  shows  a set  of  rather  small  changes  well  distributed  between 
the  four  parameters. 


In  spite  of  the  good  agreement  Just  discussed,  considerable  thought  was  given 
to  the  sources  of  error  that  might  lead  to  the  differences  between  the  experimental 
and  calculated  perturbations.  These  could  be:  a)  the  stepwise  nature  of  our  plate, 
b)  the  scale  of  our  changes  being  too  large  for  first-order  perturbation  theory , 
and  c)  errors  in  our  original  calculations  of  second  derivatives. 

With  respect  to  the  stepwise  structure  of  our  plate,  another  plate,  with 
slightly  different  design  changes,  was  built  both  in  steps  and  as  a smoothly  con- 
toured surface  interpolated  between  the  centers  of  our  sections.  The  experimentally 
determined  perturbations  showed  less  variation  between  these  two  designs  than  the 
variations  in  Table  IIIA. 


The  scale  of  our  changes,  vis-a-vis  the  range  of  first-order  perturbation 
theory,  is  the  most  likely  source  of  error.  Still,  the  errors  inherent  in  the 
calculations  of  the  second  partial  derivatives  of  the  mode  shapes  needed  to  be 
examined.  One  way  to  do  this  is  through  orthogonality  and  Hayleigh-Quotient 
tests.  To  describe  these,  let  us  define  a parameter,  K^,  as. 

s,k  ■ <*;>T  $ <*;> o6a) 

The  orthogonality  relationship  derived  in  Ref.  H requires  that 

Knk  = 0 ’ (2TA) 


and  Rayleigh's  principle  requires  that 


. i . (2SA) 

K should  be  the  modal  stiffness,  K . 
nn  n 

A program  called  CHECK  was  written  to  accept  the  data  from  the  original  uni fore 
plate  and  compute  the  checks: 


V<KnnKkk> 


1/2 


when  n/k  and 


» wheu  n=k- 


The  results  are  presented  in  Table  VA. 


TABLE  VA 


ORTHOGONALITY  & RAYLEIGH  QUOTIENT  CHECKS 


nk 


check 


12 

-.014 

13 

-.126 

23 

+.0276 

11 

0.808 

22 

0.830 

33 

0.741 

The  first  mode  does  not  seem  to  have  a stiffness  function  that  is  orthogonal 
with  that  of  the  third  mode,  and  the  three  Rayleigh  Quotients  are  all  too  small. 

It  must  be  concluded,  therefore,  that  there  are  significant  errors  in  the  second 
partial  derivatives  computed  from  the  vibration  mode  data.  This  is  understandable 
because  these  computations  involved  considerable  modeling  approximations  to  begin 
with.  From  a theoretical  standpoint,  these  checks  should  be  regarded  as  funda- 
mental requirements  for  the  modeling  of  the  mode  functions  in  the  first  place. 

Thus,  what  we  really  desire  are  values  for  the  second  partial  derivatives,  as 
close  as  possible  to  the  ones  experimentally  calculated,  that  none-the-less 
satisfy  orthogonality  and  Rayleigh-Quotient  conditions.  Fortunately,  the  minimiza- 
tion theory  presented  previously  can  be  used  to  find  these  values,  and  this  is 
presented  in  the  next  section. 

Theory  of  Minimal  Adjustments 

We  wish  to  make  the  smallest  changes  to  the  2nd  partial  derivatives  needed 
to  pass  the  orthogonality  and  Rayleigh  -Quotient  conditions. 


We  represent  the  small  changes  by  matrices, 

<o* 


(A<t>n)  to  be  added  to  the  matrices 


<«;>  i Kx.  »X,)T 

The  orthogonality  criteria  and  Rayleigh-Quotient  criteria  are  defined  by  Eqs.  (27A) 
and  (28A),  respectively.  Since  we  have  three  modes  we  have  only  six  constraints 
on  the  changes  (A4>n),  as  opposed  to  36  sections  for  each  of  three  modes.  If  we 
make  the  minimal  changes  to  the  2nd  derivatives,  however,  in  the  same  way  as  the 
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minimal  thickness  changes  were  obtained,  the  problem  becomes  uniquely  solvable  to 
a first-order  approximation.  If  we  substitute  ($n  + A^)  into  Eqs.  (27A)  and  (23A), 
the  Rayleigh-Quotient  and  orthogonality  conditions  can  be  expressed  to  first-order 

as. 


0 =/Kndxd>'  *//t(KnnVt)/V(4,n)T5(^>  * dxd>'  (as*) 

when  n=k , and 

0 ■/JPnkdxd»  * (♦;)T««p]  dxdy,  (30A) 

when  n/k,  and  AQ  is  the  area  of  the  plate.  These  criteria  can  be  multipled  by  an 
arbitrary  constant  and  added  to  the  equation  expressing  the  mean-square  change.  In 
the  case  of  two  modes , we  obtain  the  following  equation  for  zero  variance  of  the 
index  of  change. 


0 = 6//[(A*pTI  (A*’1)  + (A^)TI  (A^-y 


11 


'X2P12 


'A3P22 


] dxdy 


(31A) 


where  I is  the  identity  matrix. 

The  solution  of  tins  equation  is  found  by  setting  the  partials  of  the  integrand, 
with  respect  to  all  components  of  (A^)  and  (A$^),  equal  to  zero.  This  gives 


(A*p  = G 0]  (X),  and  (32A) 

(a^)  = g [o,  (♦Y) •(♦2)l (x)-  ;;, 


These  equations  for  (A^)  and  (A$->)  can  be  substituted  into  Eqs.  (iNA)  and  (30A)  to 
solve  for  (x).  The  three  resulting  equations  may  be  put  into  matrix  notation  to 
be  solved  for  ( X ) , given  which  A#  can  be  found  via  Eqs.  (3.?A)  and  (33A). 


s? 


Wc  wish  to  extend  this  process  to  adjust  three  modes.  Equation  (31A)  then  take; 
the  fora. 


0 *6JJ[(A*i')T  I U»f)  + (^)T  I (a*£)  + (A^)T  I (a*$) 
~X1P11  “2x2P12  *2x3P13  “X4P22  '2x5P23  'X6P33^  dxdy  * 


The  Solution;*  of  the  Euler  equation;'  in  thi.-  ease  are 


(A*p  = g [(q*),(^),(^),  o,  o,  o]  (x), 
(A^)  = G [0,  (*••),  0,  (^).(^),  0]  (X), 
(A^)  = G [0,  0, (*’'),  0 ,(*£), (*£)]  (x)  . 


( 3‘v\ ) 


( 3oA ) 

(37A) 

(3SA) 


For  simplicity,  let  us  call  the  matrices  in  the  square  brackets  of  Eqs.  (36a)-(38A), 
1.,  R,  and  S,  respectively.  Then,  the  equations  that  must  be  solved  for  \\)  may  be 
written  in  matrix  form  as 


GGQ  + RTGGR  + STGGS)dxdy]  (x). 

( 3°A ) 


As  before,  the  solution  for  (X)  is  used  in  Eqs.  (36A)-(38A)  to  generate  the  changes 
in  the  seco.  * derivatives  necessary  to  obtain  values  to  satisfy-  the  orthogonality 
and  Rayleigh- Quotient  conditions. 

A program  called  MINFUDC.E  was  written  which  performed  the  computations  indicated 
above.  Because  the  computation  was  actually  a first-order  approximation  in  Eqs.  (EdA) 
and  (30A),  two  passes  were  required  to  obtain  less  than  one  percent  error  when  the 


adjusted  second  partial  derivatives  were  put  into  the  CHECK  program.  The  result 
of  both  passes  are  presented  in  Table  VIA. 


TABLE  VIA 

CHECK  PROGRAM  RESULTS 


nk 

cneckf 1 ) 

check(2) 

12 

+ .00514 

+.000576 

13 

+.0169 

+.00285 

23 

+.00368 

+.00691 

11 

1.0306 

1.00060 

22 

1.00931 

1 . 00003 

33 

1 .0250 

1 .00024 

Finally,  the  adjusted  values  were  used  in  the  program  that  calculated  the 
perturbations  that  should  theoretically  result  from  the  design  changes  (Table  IITA,. 
Only  the  computations  via  Eq.  (ll+A)  were  made  because  these  were  presumed  to  be  the 
most  accurate.  These  are  presented  in  Table  VIIA. 


i 

TABLE  VIIA 


PERTURBATION  CALCULATIONS 


MIN FUDGE (1 ) 

MINFUDGE(2) 

Aw-|  / toi 

-.0864 

-.0822 

+ .175 

+ .174 

Aw^/ “3 

+ +.099 

+ .096 

C21 

-.2389 

-.2400 

C23 

+.0643 

+.0637 

RMS  diff. 

0.0327 

0.0328 
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It  can  be  aeon  from  Table  VITA  that  application  of  the  MTNFirpcE  program  to  the 
values  for  the  second  partial  derivatives  did  reduce  the  RMS  difference  between 
the  calculated  and  experimental  perturbations.  There  seems  to  be  little  value  in 
more  than  one  application,  however,  and  it  cannot  be  said  that  the  adjustment  of 
the  second  derivatives  made  an  outstanding  improvement  in  the  perturbation  calcula- 
tions. The  rather  large  scale  of  the  changes  attempted  in  this  design  must  be 
responsible,  therefore,  for  the  majority  of  the  errors.  It  would  be  desirable  to 
obtain  a useful  measure  of  the  range  of  application  of  this  technique.  However, 
without  a large  number  of  comparisons  between  theory  and  experiment  it  is  difficult 
to  establish  such  a range.  The  RMS  change  in  thickness,  however,  is  probably  the 
most  relevant  parameter,  and  the  value  of  28. U percent  for  the  case  presented 
here  is  probably  too  large.  Most  probably,  values  between  10  percent  and  15  per- 
cent would  be  more  in  keeping  with  the  spirit  of  first-order  perturbations. 


Finally,  some  thought  must  be  given  to  the  extension  of  this  method  to  more 
practical  structures  such  as  turbine  blades,  which  would  require  knowledge  of  the 
vectorial  mode  displacements,  and  holographic  determination  of  the  in-plane  strains 
as  well  as  the  bending  strains.  The  technique  of  heterodyne  hologram  interferom- 
etry (Ref.  9)  would  have  to  be  used  to  obtain  these  strain  values.  As  yet,  however, 
that  type  of  holographic  analysis  has  not  been  applied  to  vibration,  but,  in 
principle,  could  be  by  the  addition  of  stroboscopic  techniques.  Dividing  the  struc- 
ture into  sections  would  still  be  the  most  practical  realization  of  this  procedure, 
although  the  computation  of  the  perturbation  functions  would  be  more  complex,  of 
course.  The  simple  products  of  mode  functions  would  be  replaced  by  scalar  products 
of  vectorial  mode  functions,  and  there  would  also  be  contributions  due  to  the  in- 
plane strain  components.  The  complete  description  of  what  would  be  involved  is 
well  beyond  the  scope  of  this  paper;  however,  it  is  our  hope  that  our  elementary 
results  may  stimulate  further  work  in  this  area. 
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